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' We generalize earlier results of Fokas and Liu and find all locally analytic (f +f ^dimen- 

sional evolution equations of order n that admit an A-shock type solution with N < n + 1. 
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• d n u/dx n ) that are conditionally invariant under a given generalized (Lie-Backlund) vector 

Q , field Q(x, t, u, du/dx, . . . , d k u/dx k )d/du under the assumption that the system of ODEs 

Q = is totally nondegenerate. Every such conditionally invariant evolution system admits 
a reduction to a system of ODEs in t, thus being a nonlinear counterpart to quasi-exactly 
solvable models in quantum mechanics. 
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1 Introduction 

O 

O ', I n order to be relevant, mathematical models of natural phenomena are often required to have so- 
^J 1 ■ lutions of prescribed form (e.g. travelling wave solutions, kinks, solitons, etc.). As was first shown 
O ■ by Fokas and Liu [U [5] (see also Zhdanov [15]), a natural way to achieve this for evolution systems 
"pH . in (1+1) dimensions is to require that the system in question admit a generalized conditional sym- 
metry (GCS) such that the solutions invariant under this symmetry have the desired form. It turns 
out that the set of systems admitting GCS (such systems, or even more broadly, the systems pos- 
sessing integrable reductions, are often referred to as conditionally integrable, see e.g. [9]) is consid- 



er 

> 



^ ■ erably larger than that of systems integrable via the inverse scattering transform or direct lineariza- 
^ ■ tion, see e.g. [HE]. Notice [16] that if an (l+l)-dimensional evolution equation admits a reduction 
to a system of ODEs in the evolution parameter (time t), then it admits a GCS, and vice versa. 

This naturally leads to the following problem: how to describe all evolution systems that admit 
a given GCS? For the case of linear GCS with time-independent coefficients this problem, restated 
in the terms of the so-called invariant modules, was completely solved in the seminal paper [6] (see 
also an important earlier work [13J). For the nonlinear GCS, some partial results were obtained 
in [2], and the complete solution of the problem was obtained in [12] under certain nondegeneracy 
assumptions, see below for details; see also [6] and [12] for the survey of earlier results in the field 
and the discussion of the role played by the GCS in the search of exact solutions of nonlinear PDEs. 

However, the formulas found in [12] are not very convenient for applications. In the present 
paper we obtain alternative, easier to use, formulas for the locally analytic (l+l)-dimensional 
evolution systems that are conditionally invariant under a given analytic GCS. This is done in 
Section 3. In Section 4 we illustrate the application of our results by a number of examples. In 
particular, we completely characterize all locally analytic (l+l)-dimensional evolution equations 
of order n that admit an A^-shock type solution with N <n + l. 
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2 Preliminaries 



Consider an evolution system 

du/dt = F(x, t, u, . . . , u n ), n > 0, (1) 

for an s-component vector function u = (u 1 , . . . , u s ) T , where ui = d l u/dx l , I — 0, 1, 2, . . ., Uq = u, 
and the superscript £ T" denotes the matrix transposition. 

A smooth function of x, t, u, u\, u 2 , ... is called local (see [7]; cf. also [14] and [8]) if it depends 
on a finite number of Uj. The greatest integer m such that df/du m ^ is called the order of a 
differential function /, and we shall write that as m = ord/. If / depends only on x and t, then 
we shall assume that ord/ = 0. Unless otherwise explicitly stated, all functions considered below 
will be assumed to be local. 

A generalized vector field Q = Qd/du, where Q is an s-component local vector function is 
called [UEICS] a generalized conditional symmetry (GCS) for <Q if the system Q = is compatible 
with (JT). 

Eq-CQ) is compatible with Q = if and only if (cf. e.g. [U [15]) 

D t (Q)\ M = 0, (2) 

where D t = d/dt + ^2°1 D l (F)d /dui and D = d/dx + u i+1 d/dui are total t- and x- 
derivatives, and M. is the solution manifold of the system Q = 0. 

In what follows we shall treat Q = as a system of ODEs involving an extra parameter t, and 
assume that this system is totally nondegenerate, i.e., the systems -D J (Q) = are locally solvable 
and have maximal rank for all j = 0, 1, 2, . . ., see Ch. 2 of [8] for more details on this. Then the 
condition <^ is equivalent [8] to the following: there exist s-component local vector functions 77 - 
and an integer p such that 

s p 

A(Q) = EE^/W a )' ( 3 ) 

o=l j=o 

Finding all GCS admitted by a given system ([T]) is an extremely difficult task, whose com- 
plexity is comparable to that of the problem of finding all solutions for the system in question. 
Surprisingly enough, we succeeded [12j in completely solving the inverse problem of describing all 
systems JT]) that admit a given GCS Q = Qd/du, provided the system Q = 0, considered as a 
system of ODEs, is totally nondegenerate. 

Let us briefly review the results of [12] • Consider a system of ODEs 

Q(x,t,u,ux, . . . ,u k ) = 0, (4) 

involving i as a parameter. Here Q = (Q 1 , . . . , Q S ) T is an s-component local vector function. 
Assume that the general solution of (j3J) in implicit form can be written as 

G(x,t,u, Cl (t),...,c N (t))=0. (5) 

The number N is sometimes called the total order of (jlj). It is tacitly assumed here that G depends 
on iV arbitrary functions Cj(i), % = 1, . . . , N, in an essential way, and that det dG/du 7^ 0. 

Then using the implicit function theorem we can, at least locally, obtain the general solution 
in the explicit form as u = P(x, t, ci(t), . . . , cjv(£)), whence 

v» = d J P a (x, t, Cl (t),..., c N (t))/dx j 

Using these equations we can express c, as functions of x, tjU 1 , 
some numbers ni, . . . , , n s : 

Ci = hi(x, t,u , . . . , w ni _i, • • • , u , • • • , u ns _i), i — 1, . . . , N. 
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. ..,< 1 _ 1 ,...,m s ,...,< s _ 1 for 



Let Bi = -{dG/duy 1 dG/da,R = - {dG/du)~ l dG/dt, B { = Bi{x,t 2 h u . . . ,h N ), and 
R(x, t, u, Mi, . . . , u r ) = R(x, t, hi, . . . , hjsr), i.e., Bi and -R are obtained from Bi and R by sub- 
stituting hi for Q. Here it is understood that c« are not differentiated with respect to t while 
evaluating dG/dt. 

Let V be an open domain in the space V of variables x, t, u, u±, . . ., and let W be the set of all 
points in V satisfying the equations D^(Q) = 0, j = 0, 1,2,.. ., considered as algebraic equations. 

Theorem 1 ( [12J ) Suppose that Q(x,t,u, . . . ,Uk) = 0, considered as a system of ODEs, is an- 
alytic on V , totally nondegenerate on W , and has the same total order N on the whole ofW. 

If Q = is compatible with ([7]), i.e., Qd/du is a generalized conditional symmetry for (QJ), 
then F on V can be represented in the form 

N m s 

F = R+J2 Ut, h,..., h N )B t + E E X P , a (t, x, u, . . . , u jp jDP(Q<*), (6) 

i=l p=0 a=l 

where m andj PtCt are nonnegative integers, andQ andx pa are smooth functions of their arguments. 

By construction [12J, the system Ut = F with F given by ([6]) admits a solution of the same 
form as the general solution (jSJ) of Q — 0, that is, 

G{x,t,u,c x {t), . . .,c N {t)) = 0, 

but now the Cj(t) must satisfy 

dci/dt = Ci{t,c x , . . . ,c N ), i = l,...,N, (7) 

rather than be arbitrary functions of t. 

In other words, under the assumptions of Theorem [T] if Qd/du is a GCS for u t = F, then the 
substitution of the general solution §5§ of Q = into u t = F reduces the latter to the system of 
ODEs ©. 

3 Solving the inverse problem: which systems admit a 
given GCS? 

How to pick among the F's those of order < n, where n is a given natural number? This can 
be done using the following result. 

Theorem 2 Under the assumptions of Theorem^ on Q, suppose that 

Q a = K a -g a (x,t,u), a = l,...,s, (8) 

where g a are analytic functions of their arguments and u = (u 1 , . . . , wL_i, ■ ■ ■ , u s , . . . , ]_). 
Then the most general locally analytic F of order n>ma,x(oidR, max ord Bi, max ord hi) 

i=l,..., N j=l,...,N 

such that (QP is compatible with Q = can be locally written as 

N s n—n a _ 

F = R + £ Q{t, h,..., h N )B t + E E D m (Q a )K a>m (x, t, u, Q), (9) 

i=l a=l m=0 

where Q = {Q l ,D{Q 1 ), . . .,D n ~ ni {Q l ), . . .,Q S ,D(Q S ), . . . , D n ~ ns (Q s )) , N = & and 

K atm are arbitrary locally analytic functions of their arguments, and for n < n a K a<m = 0. 
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Proof. We have u" a = Q a + g a (x, t, u). Using this equality and proceeding inductively with usage 
of the formulas found at the previous steps, we obtain = ipf(x, t, u, Q 1 , . . . , D : '~ ni (Q 1 ), . . . , 
Q s , . . . , Di~ Us (Q s )) for j > maxn a and similar formulas for maxn a > j > n a . 

a a 

Plug the above formulas for u° into and expand the resulting expression in Taylor series 
with respect to D^(Q a ). As the order of F is n and F is locally analytic, it must be independent 
of Di(Q a ) with j > n — n a . Taking this into account, we can rewrite ([6]) in the form (Q, and the 
result follows. □ 

It is easily seen that the most general locally analytic F of order n < n = max(ord.R, 
max ord-Bj, max ordhj) such that u t = F admits a given GCS Qd/du meeting the require- 

ments of Theorem [2] also has the form with n replaced by n and K am and Q subjected to the 
extra conditions dF/duj = for j > n. 

Note that for Q = u k — g(x, t,u, . . . , u k _i) we have N = s ■ k, and Qj boils down to 

N 

F = R+J2&(t,h 1 ,...,h N )B i 

s n -k i=1 ( 10 ) 
+ E E D m (Q a )K Q , m (x, t,u,u 1} ..., UA-i, Q, D(Q), . . . , D n -\Q)). 

a=l m=0 



4 Examples 

Example 1. Let 



/3=1 j=0 



and hence iV = ^ =1 n ft . Then a general solution of Q = reads it = Eili c « where 
are arbitrary functions of £ and / \ = (//, . . . , /j S ) T are linearly independent solutions of Q = 0. 



We have HOI E] ^ = ZjZ, where 

A 1 
dfl/dx 



f! 

dff /dx 



dfl/dx 



f? 
J i 

dff/dx 



dfjf/dx 
d^^flj/dx^ 1 



fff 
df s N /dx 



d n °- L ff/dx 



n a — 1 



d n *-'f s N /dx 



n s —l 



and are obtained from Z by replacing d J ff /dx^ by w". It is easily seen that £?, = £), and 

R = J2L(Zi/z)df f (x,t)/dt. 

If • are analytic in x and £, then by Theorem [2] the most general locally analytic F of order 

n > maxordZj/Z that admits a GCS Qd/du with Q ffTTl) locally has the form 

j 

N Z- 8f (r t) N s "-"q 

^ - E -7 1' +E &(*, Zx/Z, Z N /Z)fi(x, *) + £ E D m (Q*)K a , m (x,t,u,Q), (12) 

i=l ^ CK i=l q=1 m=0 



where ^ and If a>m are arbitrary locally analytic functions of their arguments, Q = (Q 1 , D^Q 1 ), . . . , 
D n ~^(Q l ), ...,Q\ D(Q S ), . . . , D n - n °(Q s )), u = {u\ . . . , <_ x , . . . , u", . . . , < s _ x ), and K a , m = 
for n < n a . 



In the rest of this section we assume that s = 1, so for simplicity we shall write u = u, 
Q = Q, F = F. 

Example 2. Let Q = L(u), where L = Yl^ =1 (D — kj), and hi, i — 1, . . . , N, are distinct (fcj ^ kj 

if i 7^ j) nonzero constants. Then the general solution of Q = is Eili c « W exp^x), where Cj(t) 
are arbitrary functions of t, and using the results of Example 1 we obtain 

exp(-A^) 
fri = — MW, 

n (h-kj) 

where Lj = Y[j=ij^i(D ~ kj). Thus, the most general locally analytic F of order n > N — 1 such 
that Ut = F admits a GCS Qd/du with Q = can be locally written as 

N n-N 

F=J2 d{t, h,..., h N ) exp(k t x)+J2D m (Q)K m (x, t, u, Ul ,..., u N -u Q, D(Q), D n ~ N {Q)), (13) 

i=l m=0 

where Q and K m are arbitrary locally analytic functions of their arguments, and K m = if n < N . 

Example 3. Let Q = u 2 — f(u,t), where / is an arbitrary analytic function of u and t. Set 
a(z,t) = J f(z,t)dz and ip(y,z,t) = J (2a(y,t) + z)~ l l 2 dy. Then the general solution of Q = in 
implicit form reads ip(u, Ci(t),t) = x + c 2 (t), where Ci(t) are arbitrary functions of t, and we have 
h l = u\- 2a(u,t), h 2 = 4>(u, z,t)\ z=hl - x. 

Hence by Theorem [2] the most general locally analytic F or order n > 1 such that ut = F 
admits a GCS Qd/du with Q = u 2 — f(u,t) locally has the form 



F = - 

n-2 



(2a(M , t) + ftl) ^((fei>) | ^ _ Cl(t , ftl , j + _ 6(t , ftl , ft2) 



+ £ B m (Q)if m (x,t,K,m,Q,i3(Q),...,i5"- 2 (Q)), 

m=0 

where Q and are arbitrary locally analytic functions of their arguments; if n < 2, then K m = 
for all m. 

Example | As a somewhat more elaborated example, consider Q = u 2 — <p(x, t)f(ui, t), where 
ip and / are arbitrary analytic functions of their arguments, and / ^ 0. Let a(z,t) = J dzj f(z,t), 
<p(x, t) = J <p(x, t)dx, and let b(y, t) denote a solution of the equation a(z, t) — y with respect to z, 
so that a(b(z,t),t) = z. Then the general solution of Q = reads u = ci(t) + x{ x , c 2(t), t), 
where x( x > z >t) — J b{(p{x,t) + z,t)dx, q(£) are arbitrary functions of t, and we have hi = 
u - x(x,z,t)\ z=h2 , h 2 = a(ux,t) - 0{x,t). 

Hence, the most general locally analytic F of order n > 1 such that ut = F admits a GCS 
Qd/du with Q = u 2 — ip(x,t)f(ui,t) locally has the form 

^%^) +6(t, hM) + 6(t, h, h^m^l) 

dt / z=h 2 \ OZ / ^=/i 2 

n-2 

+ E ^ m (Q)^ m (x,t,n, Ml ,Q, J D(Q),..., J D"- 2 (Q)), 

m=0 

where ^ and K m are arbitrary locally analytic functions of their arguments, and if n < 2, then 
K m = for all m. 

Example 5. Assume now that Q = M(l), where M = Y[f=i(D ~~ kj ~ u )- The general solution 
of Q = has the form —v x /v, where v = Eilt 1 ^(*) exp(fcjx), and 6j(t) are arbitrary functions of 

t, cf. m. 
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Notice that —v x /v actually involves only N independent arbitrary functions of t. Indeed, 
renumbering hi if necessary, we can assume without loss of generality that h^+iif) 7^ 0. Then we 
can rewrite u = —v x /v in the form 

N 

u = -D (ln(^eyLp(k N+1 x) + y~] Cj(t) expfcx)) \ (14) 

i=l 

where Ci(t) = 6j(t)/fojv+i(0> whence 



N 

exp((k N+1 - ki)x)Mi(l) f[ (k N+1 - kj) 
i^i 

N+l 

M N+1 {1) n {h-kj) 

j=l,j=£i 



N. (15) 



,j - u) 



Here M 4 = nSW D - 

The most general locally analytic F of order n > N — 1 such that ut = F admits a GCS Qd/du 

with Q = M(l) (and hence an iV-shock type solution f|T4l) ) locally has the form 

,N+1 

N ( E {k - kj)hj exp(kjx] 

F = - (i{t, h, . . . , h N ) exp(kix) 



/N+l x 2 

(E h q exp(k q x)) (16) 

n-N 

+ E D m {Q)K m (x,t,u,u u ...,u N - u Q > D{Q),... i D n - N {Q)), 

m=0 

where Q and if m are arbitrary locally analytic functions of their arguments, and for n < N we 
have K m = for all m. This time the hi are of the form (1T51) . and for the sake of brevity we set 
h N+ i = 1. The corresponding iV-shock type solution of the evolution equation u t = F is of the 
form flT3J with that satisfy £3). 

Notice that one can readily pick among the F's (1161) those independent of x and t. They are 
of the form 

N+l 



N 7k N+1 -k N 7k N+1 ~k N ^ 22 (h kj)hj 

~r! ^ v tfciv+i-fei ' ' ' ' ' rfcjv+i-feiv-i ) /N+i „ x 2 

*— 1 "jy N 



E M (17) 

9 =1 / 



n-N 

+ E J D m (Q)^ m (M,M 1 ,...,n^_ 1 ,Q, J D(Q),..., J D"- 7V (g)), 

m=0 

where /ij = /ijexp((fcj — fcjv+i)^), and r/j and K m are arbitrary locally analytic functions of their 
arguments. 

The substitution of (114p into the equation u t = F with F of the form (TIT]) reduces this equation 
to the following system of ODEs: 

7 fcjv+i-fcjv k N+1 -k N 

aci / 62 c tv-i 



If Ci = —hi Ej==o — ^n+i)i where aj are arbitrary constants, then ffl6|) represents the most 
general F of order n such that the equation u t = F admits an iV-shock solution of the form [I] 



N+l m 

u = —D (in ^ y4j exp (kiX — t ajkj 

3=0 



where A\ are arbitrary constants. 

In other words, the class of evolution equations admitting the iV-shock solutions of the form 
(JTSJ) contains not only the equation [4] 



Ut = J2(*jD(D -u) j (u) 

j=0 

but infinitely many other equations Ut = F as well. In particular, if the order n of F is greater 
than N — 2, then Theorem [2] implies that the corresponding F's have the form 

N+l 

N m hi kj)hj 

f = E E «r(*f - fc^i) % x2 

v (?=i 7 

n-iV 

+ £ ^ m (Q)^ m (x,t,n, Ml ,...,n JV _ 1 ,g, J D(Q),..., J D n - 7V (Q)). 

771=0 

5 Conclusions and discussion 

In Theorem [2] of the present paper we completely characterized locally analytic evolution sys- 
tems ([1]) of a given order n > k — 1 that admit a generalized conditional symmetry Q(x,t, 
u, . . . , Uk)d/du with a given analytic Q of the form (JHJ). In particular, for the case of linear GCS, 
i.e., when Q is linear in Uj for all j, the right-hand sides F of such systems ([T|) are given by (fT2|) . 

The results of the present paper are somewhat easier to use in applications than the formulas 
found by us earlier in [T2]. For instance, the usage of Theorem [2] enabled us to find the explicit 
form of all locally analytic evolution equations of order n > N — 1 that admit iV-shock (and, more 
broadly, iV-shock type) solutions for a given number N, and thus generalize the results of Fokas 
and Liu [1]. 

Moreover, Theorem [2] can be employed for the classification of exactly solvable initial value 
problems along the lines of [31 |T7j and for finding evolution equations whose symmetries are com- 
patible with prescribed boundary conditions, and one can find exact solutions for the corresponding 
boundary value problems, cf. p.] and discussion in [12]. We plan to address these topics elsewhere. 
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